Quantum state preparation of a mechanical 
resonator using an optomechanical geometric phase 
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Abstract. We demonstrate that a geometric phase, generated via a sequence of four 
optomechanical interactions, can be used to increase, or generate nonlinearities in the 
£^ 1 unitary evolution of a mechanical resonator. Interactions of this form lead to new 

mechanisms for preparing mechanical squeezed states, and preparation of non-classical 
q-( states with significant Wigner negativity. 
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1. Introduction 

The geometric phase is a phase imparted on the wavefunction of a quantum state by 
driving a system around a closed path in phase space [TJ. Within quantum optics 
this phase is widely used to create logic gates for quantum computing [2J [3j H], but 
has otherwise received little attention in optomechanics. Some notable exceptions 
come from recent proposals that have considered the effect of a geometric phase 
involving mesoscopic mechanical oscillators [5J |6]. In this paper we consider an 
optomechanical system El EH HDl ttH EL2] with a time dependent optical drive that, via 
the optomechanical interactions, traverses a closed loop in phase space and thus imparts 
a geometric phase onto the mechanical element. 

It has been shown that time dependent control fields acting on an auxiliary system 
can be used to generate exotic quantum states [13]. In an optomechanical system, the 
radiation pressure force due to light in an optical resonator can be used to accelerate 
a mechanical resonator. Driving the optical resonator (being used here as the control 
field) with a suitable sequence of laser pulses can be used to manipulate the motion 
of the mechanics. In our scheme, strong mechanical non-linearity is generated with a 
sequence of four pulsed optomechanical interactions in a measurement free process. 
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During this sequence the optical field makes a closed loop in phase space and the 
mechanical oscillator obtains a phase proportional to the area enclosed within the loop 
PQ, i.e. a Berry phase. It is shown how this phase produces a nonlinear potential for 
the mechanical resonator from the linear optomechanical radiation pressure interaction, 
and in general increase existing mechanical nonlinearities. We then discuss how this 
mechanical nonlinearity can be used for quantum state preparation of the mechanical 
oscillator. 

Our full protocol takes place within a small fraction of one period of the mechanical 
oscillator, and is hence robust against rethermalisation and decoherence - similar to ref 
[8]. During the state preparation the mechanical resonator remains stationary so only 
momentum quadrature is changed. 



2. Model 

The Hamiltonian for the optomechanical interaction is given by H = hgoalaV^ib + tf), 
where go is the interaction rate, h is the reduced Planck's constant, a, (a^) and b, (Jy) 
are the annihilation (creation) operators of the optical and mechanical field respectively. 

The y/2 arises from our definition of Xm — (b + W)/y2. The equation of motion 
for a is given by a{t) = g a(t)X M — na{t) + \/2Ka in {t) where a in contains driving fields 
and vacuum noise. We now linearise about a coherent amplitude a(t) with quantum 
field operator a(t), such that a(t) = a(t) + d(t). The coherent amplitude follows the 
equation of motion a = —na(t) + y/2K0ti n (t) where a- m (t) the time dependent drive. The 
linearised field operator follows 

a(t) = go[a(t) + a(t)]X M — na(t) + \/2Ka- m (t). (1) 

The linearized approximation is made, and the a(t)Xyi term is dropped. This is 
valid even when a{t) is small as the thermal force, F t h, (e.g. in Pm = go(a + a)XM + -Pth) 
on the oscillator will still be large compared to g^aXyi- This gives an effective linearized 
Hamiltonian of the form, 

H/h = g \a(t)\ 2 X M + g \a(t)\(e ld rf + e~ w a)V2X u . (2) 

The first term generates a classical momentum imparted to the oscillator, Pm - > 
Pm + go J dt\a(t)\ 2 , where the integral is over the duration of the pulse and therefore 
proportional to the input pulse area. As this term commutes with the Hamiltonian, we 
neglect it in the following discussion and only deal with the quantum optomechanical 
interaction. If we now consider a coherent pulse input at the cavity resonance frequency, 
the unitary for the quantum interaction in equation [2] is given by, 

U{X e M ) = exp [-i X X M X M ] (3) 

where x — V%9o fdt\a(t)\ and X M = (ae~ ld + a^e l9 )/\/2. In calculating this, the 
mechanical period, Tm, is assumed to be large compared to the temporal width of 
the pulse a, so the free evolution of the mechanics may be neglected. Consider 
a pulse interacting with the optomechanical system four times; the same pulse is 
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recycled and undergoes four separate interactions. After each interaction the pulse 
is displaced in optical phase space such that the interaction Hamiltonian for each 
successive pulse is proportional to a sequence of quadrature phase operators of the 
cavity field {A M , P M , -A M , -Pm}, whh A M = (a + cfi)/y/2, P M = (a - a t )/(\/2i). For 
each interaction, it is assumed that the pulse has the same temporal profile, this is valid 
as long as the pulse width is much larger than We may use the Baker-Campbell- 
Hausdorff formula jH] to show the effective unitary is given by, 

U cS = U(-P M )U(-X M )U(P M )U(X M ) = exp[-i X 2 A M ] (4) 

where \ 2 is the area of the closed loop in the optical phase space. U e s can be seen as 
a geometric phase dependent on the position quadrature of the mechanical oscillator. 
This operator can be seen as a momentum displacement proportional to the mechanical 
position, WPmU = Pm — X 2 Xm- The momentum of the oscillator becomes correlated 
with its position, resulting in a squeezed quadrature in the final mechanical state. 

3. Experimental scheme 

Figure [1] (a) shows a schematic of one possible experiment to realise the geometric 
phase. A coherent laser pulse of temporal width a enters a large fiber cavity with 
round trip time r via a highly reflective beamsplitter. The pulse interacts with the 
mechanical oscillator via evanescent coupling from a toroidal cavity [15] with decay 
rate k, and after exiting the toroidal cavity, is displaced in optical phase space. The 
optical displacement is performed using the highly reflective beam splitter, where a phase 
controlled laser pulse is used to displace away the coherent amplitude, and displace up in 
an orthogonal quadrature Q] (b). The pulse then repeats the optomechanical interactions 
and displacement three more times to give the four pulses sequence. 

The first pulse correlates the phase quadrature of the light (Pm) with the position 
of the mechanics. The first optical displacement changes the coherent amplitude of the 
pulse from the Am (Pm) to the Pm (Am) quadrature, before the second optomechanical 
interaction. For the second interaction, Am is correlated with the mechanical position 
since it used to be Pm before the displacement. During the second interaction the 
back action of Am on the mechanical resonator effectively correlates the momentum 
of the oscillator with its position. As the position does not change over the these 
interactions interactions, correlating the momentum with the position produces a 
mechanical squeezed state. At this point the optical field is still correlated with the 
mechanical state, however, after the following two pulses the correlation is undone such 
that the final optical pulse is uncorrelated with the mechanical state, leaving the final 
state disentangled. The pulse sequence, including optomechanical interactions is shown 
in figure [TJ (b). 

To generate the geometric phase, the optomechanical system must satisfy the 
constraints, Tm 3> r > 4cr > 1/k. Setting r > 4a, ensures that each successive 
pulse decays out of the cavity prior to the next pulse entering, with only 0.01% of the 
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Figure 1. a) A schematic for an experimental protocol to realise mechanical 
nonlinearity via an optomechanical geometric phase. The required pulse sequence 
is shown entering the beam splitter, b) The evolution of the optical field in phase 
space. Blue lines show the coherent drive, with the orange lines representiong the 
optomechanical interaction. Each interaction gives a displacement in the quadrature 
orthogonal to the optical drive, with a magnitude depending on the position of the 
mechanical oscillator. 



pulse remaining. Consequently interference between successive pulses can be neglected. 
Finally we require Tm 3> t to ensure the mechanics is near motionless during the four 
pulses to generate the geometric phase. 

Figure [2] shows the effect of the four pulse sequence on the state of a mechanical 
oscillator initially prepared in the ground state, demonstrating how correlating Xu and 
Pm leads to a squeezed quadrature from a ground state Wigner function. Increasing 
the nonlinearity x benefits the protocol in two ways. Firstly it increases the effect of 
squeezing in the oscillator. Secondly it rotates the state so the squeezed quadrature 
aligns closer with the position quadrature, such that less time is required before the 
state can be verified (see section 0J, and therefore the degradation in the squeezing due 
to thermalisation will be reduced. 

4. Experimental parameters 

The previous section showed that under ideal conditions the geometric phase can be used 
to produce squeezed mechanical states. In this section we will consider experimental 
technicalities such as thermalisation of the mechanical oscillator, optical losses, and 
possible non-closing of the optical phase space loop. 

Optical losses will have detrimental effects to this protocol. The classical 
attenuation from the beam splitter will result in the phase space loop remaining unclosed 
after the four pulse sequence. This can be corrected for by changing either the amplitude 
or phase of each displacement to counteract the attenuation. We must also consider the 
amplitude-noise back action on the momentum of the mechanics. In the absence of 
vacuum noise entering each cycle, any back action on the momentum in the Xu {Pm) 
pulse will be reversed by the -Xm (— -Pm) pulse. However when vacuum noise is 
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Figure 2. a) Effect of the unitary on the Wigner function of the mechanical oscillator, 
initially prepared in the ground state. The squeezed quadrature is marked by a line in 
each graph. This quadrature is maximally squeezed at an angle tan 9 = -\/x 4 + 1 — X 2 
to the Xyi quadrature, b) The probability amplitude for the squeezed quadrature 
compared to the ground state value (shown in gray). 



introduced at the beam splitters, the amplitude noise in the — X M pulse will no longer 
perfectly cancel the amplitude noise from the Xm pulse. Unlike the attenuation of the 
classical amplitude this mechanism cannot be easily corrected for in the protocol. 

Since the beam splitters are 99:1 reflective, we expect 1 — 0.99 2 ~ 2% vacuum 
noise to be imparted onto the oscillator from each of the Xm, —Xm and Pm, —Pm pulse 
pairs. If the loss over the total cycle (beamsplitter, fiber loss, input-output coupling etc 
- moddeled as an effective beamsplitter with vacuum input) was r], then 1 — rf vacuum 
noise would be introduced to the oscillator. The square arises from the fact the pulse 
must circulate twice before it cancels the noise, e.g. the noise imparted from the Xm 
pulse will only be canceled two cycles later from the — Xm pulse. 

Even after correcting for the effect of losses, classical fluctuations in the pulse 
intensities use in the protocol could result in non-closure of the phase space loop. 
If the loop is not closed after the four pulse sequence, the effective unitary is given 
hyU = e xp[-iX M ^ jX j xti - i X 2 X M ], where X 2 = Y. 3 =i k>J xM X m ^ti}/ {2i) and 

Xi = ^9®\J XjO"iy / 7r '/2//c for a Gaussian pulse with temporal width cr, and iVj photons. 
The second in the above unitary is the geometric phase. The first term entangles the 
light with the mechanics leaving the mechanical element in an entangeled state after 
the interaction. This can be viewed as a momentum displacement on the mechanics 
that depends on the optical field, T>{— ixi oss X M loss ), where xioss is the displacement in 
the X M loss quadrature that defines the final optical state. If xioss an d 0i OS s are unknown, 
this extra phase will reduce the observed squeezing. A homodyne measurement of the 
light lost from the beamsplitter will give an estimate of xioss and 0i OS s, meaning this 
error can be accounted for retrospectively. Figure [3] shows how % loss ^ changes the 
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Figure 3. Effect of non-closure of the loop on the squeezed state for a fixed x 2 = 1. 

x loss ig 

the magnitude of non-closure in an unknown quadrature X^° ss . Once the 
variance of the squeezed state goes above one, squeezing becomes impossible as the 
amount of noise added from the non-closure of the loop is larger than the ground state 
variance. 



squeezed state. Squeezing is possible for x loss < 2.1, at this point A 2 X M > 1 meaning 
the noise added from the non-closure of the loop will degrade the observed squeezing no 
matter how squeezed the initial state is. 

Even if classical drifts in the optical displacements are corrected for, vacuum noise 
introduced by losses in the feedback loop can cause non-closure resulting in a mixed 
mechanical state. For a single pass efficiency in the fiber loop rj, the cancellation 
of noise between pairs of displacement measurements (eg. Xu and — Xm) will be 
degraded by a factor 1 — rf, leading to a loop non-closure of (expression for nonclosure 
here). For realistic inefficiencies in the range of 10%, (corisponding to 1 — 0.9 2 ~ 20% 
vacuum noise) the loop nonclosure due to non-cancellation of noise is negligible xioss = 



4<7oy 0.2(7y7r/2/K <C 1. Hence the squeezing is not significantly effected by the addition 
of vacuum noise. 

Here we consider the effect of thermalisation on the squeezed state. Thermalisation 
can have two detrimental effects. Firstly, phonon exchange with the bath during the 
four pulse sequence will render the dynamics over the pulse sequence non unitary and 
change the final mechanical state. Secondly phonons that enter during the time scale 
required for the squeezed quadrature to rotate into the measurable position quadrature 
will degrade the observed squeezing. The first of these effects can be neglected since 
the four puses can be very closely spaced with only a short delay between them. For 
example, for a mechanical oscillator with resonance frequency um = 24 kHz, and quality 
factor Q = 10 5 , the pulse duration should be a ~ 10 _8 s, such that the time for four 
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pulses (on the order of 10~ 7 s), is much smaller than the time scale for one phonon to 
enter the oscillator, 1/(IW) « 10 _5 s at IK. 

Consequently, only thermal phonons entering after the state has been prepared 
will be considered. Quantum Langevin equations were used to model the mechanical 
resonator coupled to a white noise thermal bath. The mechanical field was solved for 
the operators b(t) and tf(t) to find the variance A 2 X M = (X^) — (X M ) 2 

e r .w / „ 2 \ 



A 2 X M = —[cos(2ut)((X(,) o - (P^J + sin(2u;i) (X M P M + P M X M ) Q 

+ (*m) + (^m> - 1] + (1 - e~^)(N + 1) (5) 

where (X^} = iV + |, (P^) = (N+ §)(1 +4 X 4 ) and (X M P M + PuX M ) = -4 X 2 (iV+±) 
are the expectation values after the geometric phase has been applied. To minimise the 
initial phonon number before the four pulse sequence we envisage a cooling via pulsed 
measurement as outlined in ref [8]. In this protocol, two pulses separated by 1/4 of 
the mechanical period are used to measure the oscillator in two orthogonal quadratures, 
leading to a low entropy state. The result in j8] shows an effective thermal phonon 
number of 




1 vriV , 



where x = 4<7oy-Wp0y f I K - This gives N e g ~ 10 for a 1mm, 24kHz SiN resonator 
with Q w 10 5 and photon number Np = \a\ 2 ~ 10 6 . Using the SiN string mechanical 
oscillators considered in this paper, the effective phonon number is a achievable for 
resonance frequencies of % < 70kHz and length L < 10mm with a maximum incident 
photon flux of N p = 10 16 Hz (w 2mW at 630nm). Although this is the initial phonon 
occupation, the bath occupation remains at N ~ 10 5 . 

SiN strings present a particularly attractive mechanical oscillator, high mechanical 
quality factors of up to 7 x 10 6 have been observed, and the mechanical resonance 
frequency may be tuned via tensioning [T6l [T7] . The protocol requires the mechanical 
period to be large compared to all other characteristic time scales. From this constraint 
we will limit the following analysis to low frequency, u m = 1 — 70kHz, SiN strings. From 
|17j . the expected Q factor of a stressed SiN string of dimensions L x h x w is 

-l 



Q 



(nir) 2 Eh 2 lEh 
+ 1.0887W— - 



12SL 2 V S L 



^Bending (7) 



with E = 241 ± 4 GPa the Young's modulus of SiN, (^Bending = 17000 the quality factor 
related to bending damping mechanisms, and S = 4a;^ l L 2 psiN the tensile stress of SiN 
(with density psm) in the high stress limit. 

The optomechanical coupling rate go is calculated from evanescently coupled SiN 
string coupling rate G = 200 MHz/nm [TS] and the oscillators zero point motion: 
g = Gx = GJh/ (2mwM) with m the effective mass of the mode. The pulse width and 
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Figure 4. Progression of the variance of the mechanical oscillator over various time 
scales: a) Rethermalisation to the bath temperature, b) Oscillation of the variance 
from mechanical free evolution, c) Free evolution immediately after the four pulse 
sequence, d) Explicitly showing the squeezed region, e) Experimentally observable 
squeezing generated via a geometric phase of an oscillator initially cooled to N e g =10, 
averaged over the decay time of the cavity. The color axis gives the average variance 
in the position quadrature over one decay time of the oscillator. A value less than 0.5 
indicates squeezing below the ground state variance. The graph color axis is truncated 
at 0.5. 



optical cavity decay rate are defined by T m = 10~ 3 a = 5/k to satisfy the experimental 
requirements. 

After preparing the mechanical quantum state with the four pulse sequence, 
it may be characterised with a measurement pulse. Any measurement of the 
variance of the position quadrature will be averaged over the optical cavity decay 
time. Hence we calculate the variance measurable over this time span as A 2 X^ S = 
min T {/ 7 !" +1 / K dt'A 2 X]yi (£')}> i- e - the average variance when the oscillator is maximally 
position squeezed. Figure H] (a)- (d) shows the temporal progression of the variance over 
various time scales after the four pulse sequence. Figure 0]e) shows a plot of A 2 X^ S as 
a function of length and resonance frequency of a SiN string at IK with cross section 
157nmx3/!m. These two parameters may be chosen independently by tensioning the 
string to the desired frequency. The dashed line on the figure represents the maximum 
experimental Q of 6.9 x 10 6 p2] observed in a SiN string, with higher Q's required above 
the line. This figure shows it is experimentally feasible to achieve quantum squeezing 
for a wide range of geometries with the best squeezing of A 2 Xu ~ 10 -2 predicted for a 
3.5mm long oscillator with 20kHz resonance frequency and Q = 5 x 10 6 . For all points 
in this figure, the initial state had an effective phonon number of N e g = 10 phonons; the 
maximum intracavity photon number to achieve such cooling was 10 6 photons which is 
easily achievable. 

The X-^ appearing in the unitary is a result of the geometric phase changing a 
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linear optomechanical potential interaction into an effective quadratic potential. If the 
mechanics was instead quadratically coupled to the light field [11] [191 EH I2Q] (H = 
g ha^ aX^) , the result would be a factor X^ in the unitary - increasing the nonlinearity 
present in the Hamiltonian to forth order. In this case we may view the interaction as 
a position-cubed dependent displacement, U = exp[— ix 2 X^] = D(—ix 2 x 3 ) correlating 
the momentum of the oscillator with the cube of its position. With significant quadratic 
coupling, this provides an avenue to generate quantum states of the oscillator involving 
significant Wigner negativity see figure El 




Figure 5. Wigner functions of the mechanical state of after a geometric phase 
interaction of a quadratically coupled mechanical oscillator for values of (a) x 2 = 0, 
(b) X 2 — 0.066, (c) x 2 — 0.133, (d) x 2 = 0.2. The momentum becomes correlated 
with the cube of the position - this can be seen in the Wigner function follows a profile 
proportional to — x 3 with negativity arising in the concave sections of the curve. 



5. Conclusion 

In summary, we have demonstrated how a geometric phase in an optomecnanical system 
can be used to generate a nonlinear unitary interaction from a linear Hamiltonian, and 
in general can increase the order of non-linearity present in a Hamiltonian. This provides 
a new tool in the optomechanics toolkit allowing mechanical squeezing and quantum 
state preparation. We have shown this method to be both robust and experimentally 
feasible. 
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